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Abstract. We study the spin-down of a neutron star at- 
mosphere during the Type I X-ray burst in low mass X-ray 
binaries. Using polar cap acceleration models, we show 
that the resulting stellar "wind" torque on the burning 
shell due to the flowing charged particles (electrons, pro- 
tons and ions) from the star's polar caps may change the 
shell's angular momentum during the burst. We conclude 
that the net change in the angular momentum of the star's 
atmosphere can account for rather large frequency drifts 
observed during Type I X-ray burst. 



Key words: stars: neutron - stars: magnetic fields - X- 
rays: binaries ~ X-rays: bursts 

1. Introduction 

It is believed that thermonuclear flashes in pure helium or 
mixed hydrogen/helium layers on the surface of weakly- 
magnetic {B < 10^° G) accreting (at rate 1O~^^M0 yr~^ 
< M < 1O~*M0 yr~^) neutron stars cause the observed 
Type X-ray bursts in low mass X-ray binaries (LMXBs) 
("Hansen & Van Horn', '1975; 'Maraschi & Cavaliere, Il977 ; 
Lewin. van Paradiis. & Taam. 1995: Wooslev fc Taam . 
197fiHBildstenni99}^ . The accumulated hydrogen and he- 



lium on the surface of the neutron star by accretion ig- 
nites and burns periodically. General observed features 
of Type I X-ray bursts, such as burst fluorescences (^ 
^q39 erg), burst recurrence times (~ few hours), and 
burst durations 10 — 100 s); are successfully explained 
by spherically symmetric thermonuclear runaway mod - 
els l|Lewin. van Paradiis. k TaamL Il995t iBildstenL Il998l) . 
Recent observations from ten LMXBs with the Rossi X- 
ray Timing Explorer have revealed nearly-coherent oscil- 
lations during Type I X-ray bursts with frequencies in the 
range ~ 270 — 620 Hz (jvan der Klis, 2000: ..Strohmaveil 
l200lh . Such highly coherent oscillations are interpreted as 
rotational modulation of surface brightness asymmetries: 
burning starts locally at some points, creates a brightness 

Send offprint requests to: V. Rezania 
email: vrezania@phys.ualberta.ca 
* Present address 



asymmetry, and spreads over the entire stellar surface as 
the star rotates. This can be understood by noting that 
the timescales for accretion between bursts is much longer 
than the burst duration, so the existence of identical con- 
ditions over the whole surface of the star for the burning 
instability to start simultaneously is very unlikely. There- 
fore, these oscillations provide a direct measurement o f 
the neutron star spin frequency jStrohmaver et all Il996|) . 
Further, observing such oscillations supports the idea that 
the neutron stars in LMXBs are the progenitors of the 
millisecond radio pulsars, see iBhattacharval l)l995(l for a 
review. 

The discovery of such high coherence, large modula- 
tion amplitudes, and stable frequency oscillations while 
providing the convincing evidence that the burst oscilla- 
tion frequency is re lated to the neutron star spin frequency 
l|Strohmaveill200l[l . however, has brought many new puz- 
zles. An initial puzzle seen in the observations was that the 
oscillation frequency increases by A v ^ a few Hz during 
the burst. IStrohmaver et alJ l)l997|) firstly proposed that 
this frequency shift is due to the conservation of angular 
momentum of the decoupled burning shell from the neu- 
tron star in which the shell undergoes spin changes as it 
expands and contracts during the Type I X-ray bursts. 
Motiv ated by th i s proposal, IStrohmaver fc Markwardtl 
||1999|) and lMilleJ l|l999l l2000') modeled the observed fre- 
quency evolution and found that the oscillations are co- 
herent, as expected if rotation is the cause. But the first 
theoretical i nvestigations on understand ing the problem 
was done bv lCumming fc BildstenI l)2000|) who studied the 
rotational evolution of the neutron star atmosphere dur- 
ing a thermonuclear burst by considering one-dimensional 
vertical hydrostatic expansion. Since the nuclear burning 
time scale (> fraction of a second) is much bigger than 
the vertical sound crossing time scale through the burn- 
ing layer (~ microseconds), the outer layers of the neu- 
tron star atmosphere are always in hydrostatic balance. 
So by assuming conservation of the angular momentum 
of the shell, a hot burning layer expands hydrostatically 
by AR ^ 20 m llAvasU fc JosslllQsilHanawa fc FulimotoL 



ll984lRildst!^ri998|) , 

Av ~ Vsi2AR/R) 



and lags behind the neutron star by 
lHz(j/,/300Hz)(Ai?/20m)(10km/i?) 
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where Vg is neutron star spin frequency and R the radius. 
As the layer cools down and contracts, the rising oscilla- 
tions due to a temperature inhomogeneity, drifts upward 
as we seen by /S.v of few Hertz. By assuming that the 
burning shell rotates rigidly, they found rough agreement 
with the observed values with fractional frequency shifts 
Ai^/z/o < 0.8%. 

However, recent observations suggested that purely ra- 
dial hydrostatic expansion and angular momentum con- 
servation alone cannot account for explaining rather large 
freq uency drifts {/^v/ vn ^ 1.3%') observed in some bursts 
ICallowav ct al., 2n oitlWiinands et al.[l20nij) . As pointed 
out by IGallowav et al. il2nnih . the required expansion 
by these frequency drifts is 4-5 times large r than the 
one predicted by Gumming & Bildstcj l)200(]() . see Table 
1. In order to achieve the rather large frequency shift, 
^^^mm^^ al. ('2001*) imOTOved the calculations done by 
Gumming fc B ildstcn ( 200(1) by including the general rel- 
ativistic corrections for either slowly or rap i dly ro tating 
star. In agreement with lAbramowicz et al.l ()200lj) . they 
found that the general relativity has small effect on the 
angular momentum conservation law (~ 5 — 10%). Gom- 
paring with the data, for a rigidly rotating atmosphere, 
they obtained that the expected spin-down is a factor of 
two or three le ss than the actua l obse rved values. In an- 
other attempt. ISpitkovskv et al.l ()200 j) considered a two- 
dimensional hydrodynamics model that the burning spots 
can spread over the neutron star surface. Due to the com- 
bination of the radial expansion of the burning shell and 
rotation of the star, they proposed that horizontal hy- 
drodynamics flows may arise in the neutron star burning 
ocean during the Type I X-ray burst. By taking into ac- 
count the action of the Goriolis force due to the rapid 
rotation of the star, they showed that the horizontal flows 
may explain many features of observed bursts such as the 
very short rise time of X-ray bursts and the lack of burst 
oscillations in several bursts. Further, they argued that 
during cooling the hot ashes left by burning front, there is 
temperature gradient between equator and pole, increas- 
ing toward the pole, which drives a zonal thermal wind 
directed backward to the neutron star rotation. ^ They 
found that only if the frequency drift due to the radial 
expansion is combined with the geostrophic drift caused 
by backward zonal flows, one may expect to observe the 
rather large frequency drifts of burst oscillations in tails 
of some bursts. 

Previous models are based on the angular momen- 
tum conservation of the expanding/contracting neutron 
star atmosphere during the Type I X-ray burst. Al- 
though these models are mostly successful in explaining 
the many features of Type I bursts, the rather large fre- 



^ Actually they showed that due to the hydrostatic balance 
the material will likely ignite sooner at the equator rather than 
the pole. As a result, one would expect the burning front cools 
down at the equator firstly. 



J uency drifts observed in some bursts are still unresolved 
Strohmaver fc Bildstenl 12003^ . 

Because of no coherent pulsations seen in persistent 
emission from the majority of neutron stars in LMXBs, 
it is believed that they are weakly magnetized {B < 10® 
G). So, in previous studies the effect of magnetic field of 
the neutron star was ignored in the rotational evolution of 
burning atmosphere during the X-ray bursts. In this pa- 
per we study the change in the angular momentum of the 
burning shell during the Type I X-ray burst. This study 
is based on polar cap particle acceleration models in the 
pulsar polar cap regions. Due to the star's rotation and 
magnetic field an electric field must be induced in mag- 
netosphere, i.e. E = —{fig x r) x B/c, that has non-zero 
component alon g the magnetic field lines = E-B) (see 
for recent review lMesteill99^ . ftg is the angular velocity 
of the star, and r is a distance outside of the star. The 
parallel electric field accelerates charged particles along 
the open magnetic field lines above the polar caps up to 
ultrarelativistic energies toward the star's light cylinder 
{ricft = c). Here, we argue that, for typical magnetic field 
B ^ 10* G, the net charged particles flowing from star 
surface to infinity in the pulsar polar caps, would exert a 
stellar wind torque on the burning shell that may cause 
the angular momentum of the shell to change during the 
burst. We show that the net change in the angular momen- 
tum of the shell during the burst can account for rather 
large frequency drifts of burst oscillations. 



2. The model 

The theory of particle acceleration in pulsar magneto- 
sphere h as been studied for three de cades, after pioneering 
work by iGoldreich fc JulianI l)l969|) on pulsar electrody- 
namics. The charged particles accelerated by the parallel 
electric field (relative to the magnetic field lines) escape 
from the surface along the open field lines which extended 
beyond the star's light cylinder (of radius ric ^ 5x lO^/i^s 
cm, where Vg in Hertz) , and form a relativistic stellar wind. 
Several kind of acceleration mechanisms has been devel- 
o ped ( Meste l, 199^1, however, the one proposed bv lAronj 
(|l979lll98l|) may apply likely to X -ray burst pulsars that 
we will consider here. Arons assumed that charged parti- 
cles (electrons and ions) flow freely from the neutron star 
surface due to the thermal activities. Freely escaping par- 
ticles from the neutron star surface depends on the binding 
energy (referred as the cohesive energy for ions and as the 
work function for electrons) and the surface temperature. 
For typical magnetic field of a neutron star in LMXBs 
(■~ 10* G), the threshold temperature for thermionic emis- 
sion of electrons is T,. ~ 10"^ K and for ions is '-^ 10'^ K 
(|Luo. Shibata fc Melrosel l2000|) that are 10^ - 10^ lower 
that surface temperature of the neutron star in LMXBs. 
Further, due to the thermonuclear activities on the surface 
of the star during the X-ray burst, the particles' kinetic en- 
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ergy would increase significantly. Therefore, freely flowing 
particles is very likely from the surface of such a pulsar. 

Outward flowing charged particles above the pulsar po- 
lar caps along the open magnetic field lines causes the 
Goldreich- Julian (corotation) charge_deiisit^;_£^;y_^^^_rL • 
B/27 rc cannot be kept balanced I Harding fc Muslimovl 
Il998|) . As a result, a strong electric field develops along 
the magnetic field, above the magnetic poles due 
to the departure of total charge density from the 
Goldreich- Julian charge density, starts from zero at the 
surface and grows with distance above the surface. Even 
though pe — PG.U and therefore = at the surface, 
the curvature of the field lines causes the area of the 
open-field region to increase, so that pe increases faster 
than pa.u and then a charge deficit grows with distance 
ijHarding fc Muslimovl . Il998|) . Calculation of the parallel 
electric field is so complicated and depends on several 
mechanisms tha t work to enhance or screen t he fie l d. In a 
seri es of papers [Harding' fc M uslimov ('I998l l200l[ l200^: 
and lHarding. Muslimov fc Zhang 1,1998.) have extensively 
studied the influence of pair production, inertial frame- 
dragging effect, curvature radiation, and inverse Compton 
scattering on the E^. They showed that for the unsatu- 
rated region with altitude in range < z ^ {QsR/cY^^ — 
0.1 (for = 6007r(i/,/300 Hz) rad/s and R ~ 10^ cm) 
above the stellar surface, the parallel electric field grows 
linearly with height as E^^{9,(f>,z) = Eo{9,(j))z, while for 
the saturated region with z > (risi?/c)^/^, the parallel 
electric field that is nearly constant respect to the alti- 
tude, drops by three order of magnitudes. Here and are 
spherical polar and azimuthal angles, and z is the altitude 
in units of stellar radius. The amplitude Eq{9, (f>) depends 
on magnetic field strength B, spin frequency Vs, and ori- 
entation of the magnetic field symmetric axis relative to 
the star spin axis x- In a simple form, one can estimate 
.Ell when the charge d partic les flow fre ely from the neu- 
tron star surface (see AronsL fl^Slt lUsov fc Mehosei 119951 



iHarding fc Mushmovl l200l|) as 



Eu 2± 8.2 X 10^ 



V/m ( 



300 Hz 



5/2 



R 



10^ cm 
B 



5/2 



IQS g) (iO-4 



(1) 



Note we calculated E\\ here for a typical neutron star in 
LMXBs with nearly perpendicular field orientation rela- 
tive to the spin axis. The relative fraction of electric force 
to the gravity force above the burning front during the 
burst, ie. z ~ 10"'^, for electrons and protons are ~ 8 x 10^ 
and ~ 42, respectively. As a result, electrons, protons, and 
even ionized helium atoms can be accelerated easily along 
the magnetic field lines by parallel electric field presented 
in the polar caps of neutron stars in LMXBs. Since the 
magnetic field lines in the polar cap regions are not closed 
and extended beyond the star's light cylinder, the accel- 
erated particles in these regions may leave the neutron 



star magnetosphere. The resulting flow of particles from 
the polar caps to infinity produces a relat ivistic "stellar 
wind" , that exerts a "stellar wind torque" ijMichej Il99l|) 
on the shell, and then causes the net angular momentum of 
the shell to change during the burst. The net wind torque 
exerted on the shell due to the out flowing charged parti- 
cles from the star polar caps is estimated bv lMiche] l(l99lh 
for the Goldreich-Julian charge density to be 



Tw = /V4cf^s 



(2) 



where / — rjUQ^evA-pc is the total current of the charged 
particles flowing from the polar caps with the area ~ 
I.&ttR^/Rc that is 10% of total star's area {AttR^). Here 
riGJ - 0.36O,B/27rec ~ 2.25 x 10^^(1/^/300 Hz)(B/10** G) 
cm~^ is the Goldreich-Julian charge number density, e — 
1.6 X 10-1^ C, Rc = {GMs/nlf/^ is the corotation dis- 
tance for a star with mass Mg, and v > 0.1c is average 
speed of particles at 2 ~ 10"'^. The velocity of particle 
grows up quickly to c at z 0.1 as it is accelerated 
by 77 > 1 is a dimensionless factor (free parameter 
of the model) such that the quantity rjuoj represents the 
actual number density of particles that accelerated and 
left the star due to the parallel electric field. The value of 
?7 ~ 1 for old and isolated neutron stars (with no observed 
bursts), i.e. the space charge density is nearly Goldreich- 
Julian charge density for these stars. But in neutron stars 
in LMXBs with an ocean of accumulated hot matters, 
the value of r] must be greater than 1, due to the free 
ejection of particles from the front surface of the ocean. 
Further, during bursting the accumulated materials with 
T > 10* K on the neutron star surface , the kinetic energy 
of particles inside the ocean abruptly increases and then, 
the space charge density would increase significantly in re- 
gions with non-zero parallel electric field. As a result, one 
would expect that the number particles that may leave 
the shell would increase during the burst. 

It is necessary to note that for neutron stars in LMXBs 
with accretion rate, lO'^M© yr~i < M < lO'^M© yr~\ 
the existence of pulsar wind torque might not be consis- 
tent with accretion. One would expect the motion of the 
accreting material to oppose and shut off the radio pul- 
sar wind, and contribute to the angular momentum of the 
system. This scenario is unlikely "during the X-ray burst" 
due to the following reasons: huge thermonuclear explo- 
sion on the surface of the star causes expolsive eruption 
the accumulated materials from the surface during the 
burst with a mass bulk velocity close to the speed of light. 
As a result, one would expect the action of accretion to 
more or less shut off in this period, e.g by evacuation of 
part of the inner accretion disk during radius expansion. 
Furthermore, it has been suggested the action of accretion 
torque is balanced by the gra vitational radiation tor que 
in neutron stars in LMXBs, see lBildstenI l)l998all200^ for 
details. Therefore, the increased angular momentum due 
to accretion is lost to gravitational radiation. The latter 
mechanism provides a natural explanation for the rota- 
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tion of stars in a narrow range of frequencies, i>s 300 
Hz. Finally, the magnitude of accretion torque Tq 
M{GMsRY'^ ~ 8.6 X 10^2 dyne cm (M/lQ-^Meyr-i) 
(for Ms = IAMq and i? = 10 km) is one order of magni- 
tude smaller than that of the wind torque ^ 9.32 x 10'^^ 
dyn cm for r] - 1000 and B - 10^ G, see Eq. © and the 
discussion below Eq. 0. 

By replacing / in Eq. Q we obtain 



~ 9.32 X 10^"^ dyne cm 77 



Ms 



lAMp 



-2/3 



V300 Hz/ 



10 km/ VlO^G 



(3) 



Equation ^ gives the exerted torque on the neutron star 
surface due to the flowing charged particles from the pul- 
sar polar caps. Therefore, the burning shell of the neutron 
star will experience this torque as it expands during the 
burst. The angular momentum of the shell at a distance 
R from the center of the star in Newtonian dynamics can 
be written as 



(4) 



where constant k = {I /M i?^) shell depends on the equation 
of state and k < 2/3 (an uniform spherical shell). Here M 
and / are the mass and moment inertia of the shell. Equa- 
tion may be corrected by £ = kMR^ sin^ dflg /{I — 
2GMs/ RY^'^ for a relativistic spherical star and hy i — 
kMR^ sin^ eifls - w)/(l-2GMJR)^/^ for a rotating rel- 
ativistic star Csee lCumming et alll200ll) . Here to = Jg/R^ 
where Jj is the tot al angular momentum of the star. How- 
ever, as shown bv lCumming et all l)200lh . the relativistic 
corrections have small contributions to the final results, 
and so we ignore them. The time evolution of the angular 
momentum of the shell will be 



^(i?'l7)shcii = -T^/kM, 



(5) 



or 



(i?2fi),heii(2Ai?/i? + An/nU,n = -{T^/KM)At. (6) 

Therefore, the change in the angular velocity of the burn- 
ing shell will be 



(An/n) 



shell 



-2AR/R - (T^/nMR^n) At. 



(7) 



Here we assumed that the change in the mass of the shell 
is or der of or less th a n {AR /RY and then we neglect it 
(see ICumming et al 1 l200l|) . Equation Q represent the 
change in angular velocity of the shell during the burst. 
Now we discuss the change in angular velocity in three 
stages, before the X-ray burst or the recurrence period, 
the rising period, and after the burst. 



2.1. Before the X-ray burst 

In this period the star accumulates materials in the 
ocean that mor e or less rigidly corotates with the star 
l|Cumming et a l.. 20011. The parallel electric field inside 
the ocean is nearly neutralized due to redistribution of 
charge particles of the ocean. So in this period, the electric 
field more or less start from zero above the front surface of 
the shell and increase by height due to the charge deficit 
caused by free emission of charged particles from the front 
surface. The space charge density above the front surface 
would be bigger than but not to far from the Goldreich- 
Julian charge density and increases by time due to the 
thermonuclear activities in the oceans. So we expect that 
7? - 10 - 100. 

However, the polar cap acceleration is more or less sup- 
pressed due to strongly accreting plasma in this period. As 
a result, the pulsar wind torque would be negligible. 

Furthermore, because of nearly rigid rotation of the 
ocean with star the net wind torque during this period 
will be acted more or less on the star itself. Therefore, Eq. 
{Tj) should be corrected for this period as 



AVL 



( n 



shell 



-{T^/k'MsR^H) At, 
-4.44 X 10"^^ (-^) 



VlOO/ VI. 4M, 



Ms 



-5/3 



R 



10 km/ VSOOHz/ yiO^ G 



4/3 / B 



'0 

Atr 



104 



,(8) 



where At^cc is the burst recurrence time and k' = 2/5. As 
is clear, the net wind torque on the shell/star in this period 
and so the change in angular velocity of the shell/star is 
negligible. 

2.2. The rising time 

This period starts with abruptly thermonuclear ignition 
of the accumulated hydrogens and heliums during the re- 
currence time. The burning layers expand in les s than 
1 second up to z ^ 2 x 10"^ ( Avash k Josl Il982l: 



iHanawa fc Fufimotd.ll983lBildstelll998|) . In this period, 
as in the hydrostatic models, we assume that the rota- 
tional evolution of the expanding shell is nearly indepen- 
dent of the star itself. 

As we mentioned before, in this period due to the huge 
thermonuclear explosion on the surface of the star, the 
accumulated materials are explosively erupted out of the 
star with a mass bulk velocity close to the speed of light. 
The outgoing plasma density within period of time At ^ 1 
second, is ~ [(Apc/47ri?2)Msheii]/^pc that is 10 times (or 
more) larger than in- falling plasma density JTiacc/^pc due 
to the accretion fiow at the same period of time. Here 
"T-ace — MAt is the average mass of accreting materials 
during At. As a result, one would expect the action of 
accretion to more or less shut off during this period. 
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The parallel electric force at height z ~ 2 x 10"'^ is 80 
times bigger than gravitational force for hydrogen atoms, 
see Eq. . Therefore, the charged particles that rise with 
the expanding layers will be accelerated toward the out- 
side of light cylinder by the strong parallel electric field 
along the open magnetic field lines. Further, the sudden 
thermonuclear explosion of the accumulated materials in 
the ocean causes the thermal energy of the charged par- 
ticles to increase. Consequently, we expect that the num- 
ber density of particles that may leave the magnetosphere 
rjnQj in this period increases significantly. The resulting 
wind torque for 77 ~ 10'^ will give a net change in angular 
velocity of the shell as 

f4r) = -2\AR\/R-{T^/KMR^n) At, 

\ / shell 

(lol^) (soolfe) (io8^) (t^)'^^ 

where k — 2/3 and AR ~ 10 m is the average change 
in thickness of the burning layers during a burst. For 
77 ~ 10'^ and B ~ 10* G, the net change in angular 
velocity of the shell given by Eq. ^ is in good agree- 
ment with large frequency shifts seen in observations. We 
note that for 77 ^ 10'^, the number of particles that left 
the star AA^ ~ AirrjUQjR^ AR is much smaller than to- 
tal number of particles in the ocean N ~ M/nip, i.e. 
AN/N ~ 4 X 10"^. As expected, the change in mass of 
the shell is order of {AR/Rf - 10"^. In Table 1 we cal- 
culate the quantity 'qB for both 300 Hz and 600 Hz spin 
frequencies (for those which applicable), to produce the 
frequency shifts observed in various systems. As a result, 
the suggested magnetic fields are in range 1 — 10 x 10*" G 
for r] ~ 100 — 1000. This is consistent with other observa- 
tional evidences about the magnetic field of neutron stars 
in LMXBs. 

2.3. After the X-ray burst 

At final stage, after the thermonuclear flash, the expansion 
will be stopped by the star's gravity force, and then the 
expanded layers start to contract as they cool down. Be- 
cause of the charge redistributions during the expansion, 
the original induced electric field is neutralized inside the 
expanded shell. So the magnitude of the parallel electric 
field inside the ocean will drop to zero, and then the elec- 
trostatic acceleration will cease inside the ocean. Due to 
the lack of outward electrostatic acceleration, the ocean's 
particles are mostly accelerated downward by strong grav- 
ity. Furthermore, the action of accretion plasma is start- 
ing to recover, and as a result, the polar cap acceleration 
will shut down in this period. Although, because of new 



charge distribution, a new parallel electric field may de- 
velop above the burned front we expect that the polar 
cap particle acceleration mechanism is more or less un- 
likely after the burst. As a result, the number of particles 
that may leave the star will decrease as well, i.e. 77 < 10. 
Further as the shell cools down, it gets closer and closer 
to the star, and then its coupling to the star gets stronger 
and stronger. This can be understood by noting the mag- 
netic field near the star's surface opposes the differential 
rotation arising between shell's layers and surface of the 
star during the burst. Therefore, the net wind torque will 
be acting more and more on the whole star again rather 
than shell itself, and so we neglect it, see Sect. 12. ll 

The magnetic recoupling of the shell forces the shell 
to spin up until it achieves the star's spin frequency. The 
shell gains its angular momentum deficit (Mi?^ Af2)shGii, 
from the star, and then it reduces the angular velocity 
of the star by dVts = {M/Ms)An ~ IQ-^'^An which is 
negligible. As a result, one would expect that the oscil- 
lations' frequency will increase as the shell spins up by 
Av/vq ^ AVI/VL. In Fig. 1, we compare the corresponding 
changes suggested by previou s studies ijStrohmaver et all 
Il997t ICumming et all l200lh and the one we discussed 
here. As is clear, due to the exerted torque by the particle 
wind in the polar cap regions, the shell's spin frequency, 
and then the oscill ations' frequency decrease more than in 
the calculation bv lCumming et all l|200l|) during the ris- 
ing time. Therefore, we expect to observe larger change in 
oscillations' frequency during the burst tail. 

3. Discussion 

Among the ^ 50 known Galactic LMXBs, the highly co- 
herent burst oscillations with large modulation amplitudes 
and stable frequencies in range ~ 270 — 620 Hz are seen 
in ten LMXBs during Type I X-ray bursts. These oscilla- 
tions are most commonly seen during tails of bursts, when 
the burning is thought to have spread over the whole sur- 
face and obviously asymmetry is no longer present, are 
not observed in all Type I X-ray bursts from the same 
source. While it is believed that the oscillation frequency 
is closely related to the neutron star spin frequency due 
to the observation of the kHz quasi-periodic oscillations in 
the persistent emission ()van der Klij . l200fil) . there is still 
an ambiguity that whether the oscillation frequency is the 
spin frequency or twice the spin frequency. Further, as seen 
in observations, the oscillation frequency increases by Av 
a few Hz d uring the burst. Th i s freq uency shift is firstly ex- 
plained bv lStrohmaver et al] l|l997|) that the burning shell 
decouples from the star, and undergoes spin changes due 
to the conservation of angular momentum of the shell as 
it expands and contracts during the Type I X-ray bursts. 
Further observations and studies suggested that purely ra- 
dial hydrostatic expansion and angular momentum conser- 
vation alone cannot explain rather large frequency drifts 
[Av/v ~ 1.3%) observed in some bursts ()Cumming et al.l . 
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l200l[lrTallowav et allEoOltlWiinands et alll20nij) . For re- 
cent revi ew on the thermo nuclear bursts and their prop- 
erties see IStrohmaver iST^ ildstcn ( 2003). 

In this paper, we addressed the latter problem by 
studying the evolution of angular momentum of the burn- 
ing shell during the Type I X-ray bursts in LMXBs. Based 
on p article accelera tion models near a pulsar polar cap re- 
gion l)Mesteifl998l) . we studied the change in the angular 
momentum of the burning shell during the Type I X-ray 
burst. The net charged particles that accelerated by par- 
allel electric field in star's polar caps, flow from the star 
surface to infinity through open magnetic field lines, would 
exert a torque on the burning shell (called wind torque) 
and cause the angular momentum of the shell changes dur- 
ing the burst. We introduced a dimensionless factor 77 that 
the quantity rjncj represent the number density of the 
ejected particles from the burning layers above the pulsar 
polar caps. We showed that for -q ~ 100 — 1000 (during 
X-ray bursts) and a typical magnetic field B ~ 1 — 10 x 10^ 
G in the pulsar polar caps, the rather large observed fre- 
quency drifts of burst oscillations can be explained by the 
resulting wind torque exerted on the burning shell. In Ta- 
ble 1, we obtained the values for tjB for both 300 Hz and 
600 Hz spin frequencies that causes the corresponding fre- 
quency drift observed in each burst. As is clear, the result- 
ing magnetic fields' strength is in a good agreement with 
other observations from the neutron stars in LMXBs, such 
as the lack of coherence pulsations in persistent emission. 

Above the neutron star's polar caps particles flow out- 
ward along the open magnetic field lines and a steady 
charge density cannot be maintained at the Goldreich- 
Julian density everywhere. Consequently, a strong elec- 
tric field develops along the magnetic field, extracts the 
'primary' charged particles with density n ^ uqj (for 
a star with no burst activity, see Sect. [JJ from the sur- 
face and accelerate them to high Lorentz factor (~ 10^). 
According to the polar cap models, the predicted radi- 
ation luminosity for a typical millisecond pulsar would 
be £ ~ 1033(z/,/300 Hz)5/28(S/108 G)-i/7 ergs/s that 
is comparable to 7-ray luminosities. So, one may expect 
to observe 7-rays in millisecond pulsars that h as been de- 
tecte d only in PSR J0218-I-4232 by EGRET llKui et. all 
l20onh . However, polar cap particle acceleration is limited 
by various energy loss mechanisms such as curvature radi- 
ation, resonant and nonresonant inverse Compton scatter- 
ing, and pair production. The accelerated particles along 
the magnetic field lines suffer energy loss and emit photons 
through the curvature radiation. Further, these particles 
scattered by soft thermal X-ray photons emitted from the 
hot surface and loss energy via inverse Compton scatter- 
ing. In addition, the high energy photons resulted from 
the curvature radiation and/or inverse Compton scatter- 
ing produce a cascade of 'secondary' pair particles 
with a number density n^± ^ lO^n above polar cap re- 
gion that change the charge distri butions and screen the 
electric field (see for more detail iHarding fc Muslimovl 



. In millisecond pulsars (with no burst activities) 
the po lar cap particle acceleration m echanism was stud- 
ied bv lLuo. Shibata fc Melros^ l|2000l) . They showed that 
in these pulsars this mechanism, which is relatively effi- 
cient in accelerating charged particles in the star's polar 
caps, mostly suffers energy loss from the radiation of ac- 
celerated particles along the field lines ( curvature radi- 
ation). So the maximum energy attainable by the parti- 
cle is limited by energy loss through the radiation reac- 
tion. They found that the maximum Lorentz factor that 
might be achieved by an individual particle at z ^ 0.1 
is r = {GneoE^^p^/eY^^ ^ 10^, however, the Lorentz fac- 
tor for the bulk pair plasma is much smaller ~ 100. The 
latter may explain the lack of detection of high energy 7- 
rays in these stars. Here is given by Eq. |^ and pc ^ 
(4/3)(ci?/27r:/^)i/2 _ 5 X 10^^(300 Uz/iysY^^iR/lO^ cm) m 
is the curvature radius of the field lines. 

For neutron stars in LMXBs with intermediate accre- 
tion rate 10""Mq yr"i < M < IO'^Mq yr"\ how- 
ever, the polar cap acceleration mechanism mostly suf- 
fers from the action of accretion plasma. Because of the 
large plasma density (wacc/^pc) of the accretion flow to- 
ward star's polar caps, the parallel electric field in the po- 
lar cap regions would be shorted out and the polar cap 
acceleration mechanism would be suppressed. We note 
that, this scenario is not likely during X-ray bursts due 
to the explosively erupting materials from the surface 
and possible evacuation of part of inner disk. In this pe- 
riod with At ~ 1 second, the eruptive plasma density 
is 10 times larger that the accreting plasma density, ie. 

- [(Apc/47ri?2)M^j^^ll]/^p^ ^ 10 maccMpc 
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Table 1: In this table we obtain the lowest radius expan- 
sion AR (based on hydrostatic expansion models) and the 
scaled magnetic field -qB (based on the polar cap parti- 
cle acceleration models) for some X-ray bursts in LMXBs. 
We calculate these quantities for two chosen values of the 
star spin frequency, i^a ~ 300 Hz and Vs ~ 600 Hz, such 
that the corresponding frequency shift is comparable with 
observations, vq is the oscillation frequency and Av its cor- 
responding shift seen during the burst. The value of rjB is 
based on the assumed value of AR ~ 20 m for both spin 
frequencies, see Eq. 

Table 1 



Object 


Time 


fo 


Ai^/uo 




(UT) 


(Hz) 


(10-3) 


4U 1636-54 


1996 Dec 28 (22:39:22)'''' 


580.5 


~ 2-4 




1996 Dec 29 (23:26:46)^'^ 


581.5 


~ 3 




1996 Dec 31 (17:36:52)^'^ 

1997 Jul 26 (14:04:19)* 
1997 Jul 30 (12:11:58)'' 

1996 Feb 16 (10:00:45)'*''' 

1997 Sep 9 (06:42:56)'* 


581. 


~ 3 


4U 1702-43 


329.85 ± 0.1 


7.7 ± 0.3 




330.55 ± 0.02 


4.8 ± 0.3 


4U 1728-34 


364.23 ± 0.05 


6.6 ± 0.1 




364.10 ± 0.05 


5.9 ± 0.2 


Aql X-1 


1997 Mar 1 (23:27:39)®''' 


549.76 ± 0.04 


4.3 


4U 1916-053 


1998 Aug 1 (18:23:45)* 


269-272 


13.0 


MXB 1658-298 


1999 Apr 14 (11:44:52)^ 


567 


9.0 
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IStrohmaver et alJlll996l): ('6llZhang et 3 .1.71998^: (7)'Fox et _al' 

QML); fSl .Galloway et all (|200J) : Cg HWiinands et al. L2 ^ 
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Fig. 1. In this figure we schematically compare the change in angular velocity of the expanding shell and then 
oscillations' frequency during the type I X-ray burst suggsted by ICumming et alJ l|200l|) (dotted line) and the one 
proposed here (solid line). Due to the exerted torque on the shell by the particle wind from the magnetic po l ar ca p 
regions during the burst, we expect a larger depth in Afi rather than the model discussed bv lCumming et all |200l|) . 
As a result, the larger drift would be expected in the burst tail. The solid and dotted curve after the burst are drawn 
schematically to show star-shell recoupling in this period. 



